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We have used the Australia Telescope Compact Array to seek the integrated radio 
flux from all the pulsars in the core of the globular cluster 47 Tucanae. We have 
detected an extended region of radio emission and have calibrated its flux against the 
flux distribution of the known pulsars in the cluster. We find the total 20-cm radio 
flux from the cluster's pulsars to be S* = 2.0 ± O.SmJy. This implies the lower limit 
to the radio luminosity distribution to be ™*"Li4oo = 0.4mJy kpc and the size of the 
observable pulsar population to be iV< 30. 



1 INTRODUCTION 



Pulsars fall into two quite distinct categories: young objects 
with relatively low spin rate — standard pulsars; and much 
older objects spinning several hundred times per second — 
the so-called millisecond pulsars. The majority of known pul- 
sars fall into the former category and are distributed across 
the Galactic disk as expected for objects that have origi- 
nated fairly recently (10* - 10^ years) from massive stellar 
supernovae. The millisecond pulsars have quite a different 
distribution, but which can also be understood by looking at 
their progenitors. Srinivasan & van den Heuvel (1982) and 
Alpar et al. (1982) proposed that the millisecond pulsars 
have been spun up to high rotation rates during an accre- 
tion phase as a member of low mass X-ray binary (LMXB) 
systems. The tendency for millisecond pulsars to be found in 
globular clusters is seen as a natural consequence of the high 
incidence of LMXB systems in such clusters. Of the galactic 
globular clsuters, 47 Tucanae dominates as the host of at 
least 22 pulsars (Manchester, 2000). This number is large 
enough that it presents the possibility of certain population 
studies, some of which have been addressed by Camilo et al. 
(2000) and Friere et al. (2001). 

In this paper we discuss the number of potentially ob- 
servable pulsars in 47 Tucanae. We use the term "poten- 
tially observable" (sometimes simply "observable") for pul- 
sars whose emission beam passes over the Earth and so could 
be detected with a sufficiently sensitive telescope. The prob- 
ability that any pulsar is potentially observable depends on 
the angular width of its emission beam, and the fraction of 
all pulsars that are observable relates to some population 
average (latitudinal) beam width. 

Camilo et al. have estimated the total number of poten- 
tially observable pulsars in 47 Tuc to be ~200. They base 
their estimate on the supposition that the cluster pulsars are 



distributed in luminosity according to the law observed for 
the standard pulsars in the Galactic disk. That is, their dif- 
ferential luminosity distribution is assumed to be consistent 
with 



dlogiV = -dlogL 



(1) 



where A'' is the number of pulsars and L their luminosity. 
Camilo et al. also assumed that the pulsars extend in lu- 
minosity down to the lower limit observed in the Galac- 
tic pulsar population, measured at 400 MHz, of '"'"L400 ~ 
1 mjy kpc^. 

The determination of the luminosity distribution for a 
population of pulsars is a difficult process (see for example 
Lyne, Manchester & Taylor, 1985) for the following reasons. 
Firstly, the distance to each pulsar in the sample must be 
known. This is usually inferred from the observed disper- 
sion of the pulsar signal and estimates of the free electron 
density along the interstellar line of sight and is often given 
with large uncertainty. Secondly, any sample is incomplete 
because of observational selection effects. This is particularly 
relevant for the fainter pulsars in the population and intro- 
duces extra uncertainty in the measured luminosity function 
at the low- luminosity end of the distribution. 

A measurement of the total radio flux from the core 
region of 47 Tucanae would, if the form of the pulsar flux 
distribution were known, provide an estimate of the lower 
limit to the distribution and of the number of pulsars in the 
cluster. Specifically, if the cluster contains the estimated 200 
observable pulsars and if the flux distribution as assumed by 
Camilo et al. is followed, the total radio flux from the pulsars 
in the cluster core should be ~4 mJy — double that of those 
already detected. Here we report observations designed to be 
sensitive to such a distribution of 20 cm emission from the 
core. The observations reported here were made with the 
Australia Telescope Compact Array (ATCA) in low resolu- 
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tion (short baseline) configurations and are complementary 
to earlier higher resolution ATCA observations reported by 
McConnell & Abies, 2000 (referred to hereafter as MA2000) 
and McConnell, Deacon & Abies, 2001 (MDA2001). 



2 THE PULSAR LUMINOSITY 
DISTRIBUTION 

The observed radio luminosity of a pulsar is 



L = sD^ 



(2) 



where s is its radio flux and D its distance from the observer. 
L is usually expressed in the units mjy kpc^. The distance 
to 47 Tuc is 4.5 ± 0.3 kpc (Zoccah ot al., 2001). In equation 
(2), the distance D to all pulsars in 47 Tuc is effectively the 
same. Given the angular spread on the sky of the 47 Tuc 
pulsars A6 ^ 2'.5, and assuming the distribution to bo sym- 
metrical, the spread in pulsar distance AD DA9 ~ 3pc 
is negligible. In the following we describe the distribution of 
radio fluxes, and in conclusion give the equivalent parame- 
ters of the luminosity distribution for comparison with the 
general pulsar population. 

As noted in the introduction, the population of stan- 
dard pulsars in the Galactic disc has a luminosity distribu- 
tion consistent with equation 1. It is useful to write explic- 
itly the underlying frequency distribution and to determine 
whether the 47 Tuc pulsars with known fluxes also belong 
to a population consistent with that relation. We write the 
frequency distribution of pulsars, n{s), with radio flux s as 



n{s)ds ■ 



2iV SminS 





S <C Smin 



where 



f 



n{s)ds = N 



(3) 



(4) 



is the total size of the pulsar population. The cumulative 

frequency distribution A'^(> s), ic. the number of pulsars 
with flux density greater than s, is then 



N{> s) 



f 

J s 



n{x)dx ■■ 



NSr, 



(5) 



for s > Smin ■ After taking the logarithm we get 

log N{> s) = log NSmin " log S (6) 

which, in differential form, corresponds to equation 1. For 
such a distribution, the integrated flux of all pulsars with 
flux above some limit s is 



5'(> S) = NSmin ^ i 
fc=l 



(7) 



where m — Ns,nm/s. 

To verify that equation 6 describes the sample of pulsars 
already detected in 47 Tuc, we have analysed their fluxes re- 
ported by Camilo et al. (2000) . We give their cumulative flux 
distribution A'^(> s) in Figure 1, which shows that pulsars 
with s> 0.07 m,Jy follow the assumed law. Fainter pulsars 
appear to be less numerous, as might be expected either 
because of approaching the sensitivity limit of the detec- 
tion equipment or because the intrinsic lower flux limit is 
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Figure 1. Tlic cumulative distibution of 47 Tuc pulsars with 
respect to 20 cm radio flux taken from Camilo et al. (2000). The 
line is a fit to the distribution for pulsars with s ^ 0.07mjy, 
constrained to have slope 6 = —1: logAr(> s) = —0.18 — logs. 



being approached. We have fitted the cumulative distribu- 
tion for s ^ 0.07 mJy (9 pulsars) to an expression of the 
form log A'^(> s) — a + b log s and find that a = —0.1 ± 0.2, 
b = —0.9 ±0.2, consistent with equation 6. Constraining the 
distribution to follow equation 6 exactly {b = —1) we find 



log iV(> s) = -0.18 ± 0.07 - log s 
This indicates that the quantity Nsj, 



(8) 

0.66±0.11mJy. 



3 OBSERVATIONS AND DATA ANALYSIS 

The aim of the observations reported here was to detect 
radio emission from the core of 47 Tucanae. The likely dis- 
tribution of pulsars in the cluster has already been indicated 
by the positions of the brighter members as determined from 
pulse timing analysis (Preire et al., 2001). We sought to de- 
tect the integrated emission of a possibly larger population 
of pulsars sharing this spatial distribution. All pulsars lie 
within 3.2 Vc of the cluster centre, where the core radius 
rc = 23'.'1 ± 1'.'7 (HowcU, Guhathakurta & Gilliland, 2000). 
We therefore chose to use the Australia Telescope Compact 
Array in the 0.75 and 0.375 configurations to give sensitivity 
to brightness distributions of order ~50" in size. The obser- 
vations were made over a bandwidth of 128 MHz centred at 
1.408 GHz (A ~ 20 cm). 

The MIRIAD data reduction package was used to cali- 
brate the data and form images. The antenna gains were 



Table 1. Summary of the ATCA observations of 47 Tucanae. 



Date 


Time 


Integration 


Array 


(UT) 


(UT) 


(hours) 




2000 Aug 05 


14:16 


9.0 


0.375 


2000 Dec 24 


00:10 


11.2 


0.75C 


2000 Dec 25 


23:13 


18.7 


0.75C 


2000 Dec 27 


11:59 


8.9 


0.75C 



calibrated using brief observations of the source B2353-696 
(5*1.4 = 1.05Jy). The flux scale was referred to the ATCA 
primary flux calibrator B1934-638 (5i.4 = 14.9 Jy). Some 
improvement to the antenna gain calibration was achieved 
by self calibration of the final image. 

The images we present here were derived from the ob- 
servations listed in Table 1. Care was taken to image a large 
area (137' x 137') to cover the annular region of sky falling 
in the first sidelobe of the 22m antenna beam pattern which 
lies about 60' from the centre of the primary lobe. A number 
of sources with flux density of up to ~6 mjy were detected 
in the first sidelobe. 



4 IMAGE ANALYSIS 

The image is dominated by the many background sources 
in the field, the brightest having a flux of ~50 m,Jy and 
lying 6'.5 from the image centre. The effect of these sources 
(including their possible variability), combined with unmod- 
elled antenna gain variations, is to increase the image noise 
from the 20/xJy/beam expected from thermal (receiver and 
sky) fluctuations to the measured value of ct = 60/iJy/beam. 
Nonetheless, a diffuse area of 20-cm emission is visible in the 
cluster core region of the image shown in Figure 2. 

Before measuring the flux of that emission and attribut- 
ing it to the cluster pulsars, we must consider a number 
of effects associated with flux calibration, the sensitivity of 
the array to extended emission and pulsar variability, all of 
which may influence the measurement we wish to make, and 
its interpretation. 



4.1 Sensitivity to extended emission 

The expected spatial extent of pulsar radio emission is in- 
dicated by the radial distribution of the 15 pulsars with 

locations derived from pulse timing analysis (Friere et al, 
2001). These pulsars, to an approximation sufficient for our 
purposes, are uniformly distributed in radial distance from 
the cluster centre out to 3.2rc. Thus half the pulsar radio 
flux is expected from the inner quarter of a circle of ra- 
dius ^ 75". To test the sensitivity of the ATCA to the pul- 
sar fiux we constructed a simulated visibility dataset using 
the same antenna locations, observation hour angles and re- 
ceiver noise characteristics as those of the real observations. 
To these data we added contributions corresponding to: (i) 
a circular source with total flux 1 mJy and a gaussian pro- 
file of full width to half maximum of 75"; (ii) a set of 17 
point sources whose fluxes (each <100 nJy, total 946 fJ,3y) 
were drawn from the distribution described in section 2 and 
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Figure 2. 20-cm image of 47 Tucanae. The circle is centered on 
the cluster centre and has radius 4 x rc where Tc = 23'.' 1. The 
size of the restoring beam is shown at lower left. Image brightness 
is shown as contours at -3,3,5 & 8 X the rms fluctuation in the 
image of <t = 60/^Jy. 



whose positions were drawn from the uniform radial distri- 
bution mentioned above. In the images derived from both 
these trial datasets the measured flux of the added com- 
ponent differed from the expected value by less than the 
standard error expected from the simulated thermal noise. 
We conclude that the baseline lengths used in the obser- 
vations adequately sample the spatial scales present in the 
likely pulsar emission. 



4.2 Flux calibration 

We aim to compare the radio fiux of pulsars mesisured with 
the Parkes radiotelescope (Camilo et al., 2000) against the 
flux of radio sources in images formed with the ATCA. 
Therefore it is necessary to ensure consistency between the 
flux scales used in the two sets of measurements. In the 
case of the Parkes observations the pulsar time series sig- 
nal strength is compared to the amplitude of noise fluctu- 
ations arising from components of the receiving equipment 
and from the sky — the so-called system equivalent flux den- 
sity. Camilo et al. estimate a systematic uncertainty of or- 
der 25% in their published flux scale. The amplitudes of the 
visibilities measured with the ATCA were calibrated using 
observations of the source B1934-638. No specific attempt 
has been made to tie together the flux scales of the two 
instruments (F. Camilo, private communication). 

To make a direct comparison we use flux measurements 
of the five strongest pulsars (C, D, E, F and J) with both 
instruments: the Parkes ineasm-ements are those of Camilo 
et al. (2001); the ATCA measurements are from the earlier 
high resolution observations (MDA2001) which, unlike the 
observations described here, used the 6-km conflgurations 
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of the ATCA to achieve its highest possible resolution. The 
measurement results appear in Table 2. We use the symbol S 
for integrated fluxes, distinguishing between the bright (five 
brightest) and faint (sixth brightest and fainter) pulsars with 
the leading superscripts b and /; hence ''Sp is the sum of 
fluxes for pulsars C, D, E, F and J measured at Parkes, 
and ^Sa is the ATCA measurement of flux from all fainter 
pulsars. We calculate the total fluxes: 

''S'p = 1.48 ±0.09mJy (9) 

''S'a = 1.20 ±0.16mJy (10) 

To validly compare fluxes measured in the ATCA images 
with the flux distribution quoted by Camilo et al. we scale 

the ATCA data by the factor ''Sp/''Sa = 1.23 ± 0.18. This 
factor is consistent with the Parkes flux scale uncertainty of 
~25% estimated by Camilo et al. 

4.3 Pulsar variability 

A remarkable feature of the 47 Tuc pulsars is their strong 
scintillation. Camilo et al. (2000) report that pulsar fluxes 
measured over a series of 17. 5- minute observations may vary 
by factors of up to ~30. In MA2000 we reported variability 
of the imaged pulsars over a series of 12-hour integrations. 
We have extended that analysis using the additional ATCA 
observations described in MDA2001 and find that the mea- 
sured fluxes of an individual pulsar from a set of ~ 12-hour 
integrations has a standard deviation approximately equal 
to the mean flux of that pulsar. Note that for diffractive 
scintillation, the expected exponential distribution of fluxes 
(Rickett, f977) has standard deviation equal to the mean, 
consistent with this observation. Thus, even over the 48-liour 
integration of Figure 1, scintillation of the brightest pulsars 
in the cluster could significantly infiuence the integrated fiux 
from the core and take us to a misleading conclusion. 

It is relevant to consider the variability of total flux of 
the cluster's pulsar population. Consider a set of variable 
sources each with mean flux Si and rms fluctuations in flux 
of (5si ~ / X Si after integration over time t, where / is 
fixed for the whole population. Then if the sources have a 
luminosity function as given by equation (1), it follows that 
the rms fluctuation of the total flux is 

AS ~ 1.3/ (11) 

where Smax is the flux of the brightest member of the 
set. In the case of the 47 Tuc pulsars / ~ 1 for 12-hour 
integrations, so we would expect 12-hour measurements of 
the total flux to fluctuate with AS* ~ 0.7mJy. This quantity 
is dominated by scintillations of the brightest pulsars in the 
cluster. Independent measurements of their individual fluxes 
over the integration period considered would improve the 
estimate of the mean total flux, provided the measurement 
error is less than the expected rms fluctuation over the same 
period. 

The imago presented in Figure 1 was observed with five 
elements of the ATCA arranged to give baselines of up to 
750m, chosen to give best brightness sensitivity to extended 
emission, but which naturally lacks the angular resolution 
to measure fluxes of individual pulsars. However, the ATCA 
has a sixth element fixed permanently 3km from the rail 
track supporting the moveable five, so that each observation 



included visibility measurements from five baselines with 
lengths in the ranges 4.3 - 5.0 km (0.75C array) or 5.5 - 
6.0 km (0.375 array). The image derived from all baselines, 
inspite of the poor sampling of the uv-plane, has a synthe- 
sised beam of 4'.' 1x3'.' 5 and so is sensitive to point sources. 
The image has noise fiuctuations of cr = 40/jJy/beam and all 
of the brightest five pulsars are visible at > 3cr (Table 2, col- 
umn 4) . Thus flux measurement of the five brighter pulsars 
introduces significantly less uncertainty (~ y/5 x 0.04 mjy) 
than expected from the scintillations. 



4.4 Image quality in the presence of variable 

sources 

The fidelity of radio aperture synthesis images depends on 
the constancy of the region being imaged. In reality vari- 
able sources exist and some degradation of any 20cm image 
from the ATCA can be expected. For these observations we 
expect the pulsars themselves to vary. The diffraction re- 
sponses around a variable source will differ from the point 
spread function of the observation and so will not be mod- 
elled correctly during deconvolution. In aperture synthesis 
the point spread function has zero mean, so that the effect 
on the image some distance from the varying source will 
be as often positive as negative. The effect of a number of 
independently varying sources in the field is to increase the 
background fluctuations in the image. The images presented 
here have rms brightness fluctuations a — 60/iJy/beam, 
compared with the value a = 20/iJy/beam expected from 
thermal considerations alone. We can attribute at least part 
of this extra noise to variations in the background sources. 

To gauge the size of fluctuations induced by source vari- 
ability we have artificially added a 0.5mJy source to the ob- 
served data over a 11-hour section of the 48-hour integration 
and determined the effect on the fiux measured in the core 
of the cluster. With the artificial source placed in the core, 
the measured flux is increased by the expected time aver- 
aged flux of ~0.12 mJy. If the test source is placed outside 
the core region, the infiuence is position dependent, but the 
change to the measured core flux is < 25/xJy. 



4.5 Analysis method 

To determine the best possible estimate of the total pulsar 
flux, on a scale that allows direct comparison with the flux 

distribution analysed in section 2, we have analysed the data 
in the steps described below. 

• Measure the fluxes for all pulsars above a 3a detection 
limit in the high resolution image described in section 4.3. 

• Subtract from the low resolution visibility data, a set 
of point source responses with the measured fluxes at the 
corresponding pulsar locations. 

• Subtract from the low resolution visibility data all 
source components detected in the high resolution image 
outside the circular region centred on the core with radius 
4rc. 

• Subtract an additional source (S176 in MDA2001) 
which lies within 4rc of the cluster centre. This source was 
discussed in MDA2001. It lacks variability, has a flat spec- 
trum and is not suspected of being a pulsar. 
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Table 2. Flux measurements of the pulsars in 47 Tucanae. (2): 
fluxes measured at Parkes (Camilo et al., 2001). (3): fluxes mea- 
sured from the ATCA image reported in MDA2001. (4): fluxes in 
the image formed from the long baseline data as part of the obser- 
vations listed in Tabic 1. Uncertainties in the last digits quoted in 
(2), (3) appear in parentheses. The uncertainties of values in (4) 
are ~ 0.04 mjy. 



Pulsar 


Spks 


Satca 


Si ATCA 


name 


(mJy) 


(mJy) 


(mJy) 


(1) 


(2) 


(3) 


(4) 


C 


0.36(4) 


0.24(6) 


0.13 


D 


0.22(3) 


0.37(10) 


0.22 


E 


0.21(3) 


0.08(2) 


0.14 


F 


0.15(2) 


0.18(5) 


0.14 


G 


0.05(2) 






H 


0.09(2) 






I 


0.09(1) 






J 


0.54(6) 


0.33(9) 


0.16 


L 


0.04(1) 






M 


0.07(2) 






N 


0.03(1) 






O 


0.10(1) 






Q 


0.05(2) 






u 


0.06(1) 







• Derive the low resolution image from the modified vis- 
ibility dataset, and measure the flux of all emission within 
Arc of the cluster centre. 

• Correct this figure for the flux scale diferences between 
the ATCA and Parkes measurements (see section 4.2). 

• Add the mean flux,as reported by Camilo et al. (2000), 
of each pulsar detected and subtracted in the first two anal- 
ysis steps above. 



5 RESULTS 

Figure 3 shows the central part of 47 Tucanae at 20 cm be- 
fore (left) and after (right) the pulsar and source removal 
described in section 4.5. Significant emission is visible from 
only a small part of the Arc circle around the cluster centre 
and the total radio fiux is much less than the ~4 mJy ex- 
pected if the predicted ~200 pulsars exist in the cluster. The 
total flux measured in the residual image inside the Arc circle 
is ^ Sa = 0.42±0.22m,Jy. The uncertainty is calculated from 
the observed brightness noise in the image (0.06 mjy/beam), 
the error in determination of the bright pulsar flux, and an 
estimate of the likely difference in flux of a 48-hour inte- 
gration from the long term mean flux of the weaker, un- 
subtractcd pulsars (see section 4.3 and equation (11)). Con- 
verting the measurement to the Parkes flux scale using the 
factor derived in section 4.2, we have ^ Sp = 0.52±0.28mJy. 
Thus the total 20-cm radio flux is 

Sp =^ Sp +^ Sp = 2.00 ± 0.29mJy (12) 

In Figure 4 we relate this result to the cumulative 
flux distribution. The vertical line indicates the likely lower 
bound to the flux of pulsars in the cluster. It is unlikely that 
the flux distribution extends below Smin — 0.02 mJy, corre- 
sponding to an upper bound for the pulsar population size 
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Figure 3. Image of the core region of 47 Tucanae. The left panel 
shows the radio emission from all pulsars and background sources. 
On the right is the residual image after subtraction of components 
outside the Arc circle, pulsars C,D,E,F and J, and source S176 
(see text). Image coordinates arc indicated as offsets from the 
cluster centre: aj2000 = 00 : 24 : 05.9, 5j2000 = -72 ; 04 : 51.1. 
Image brightness is shown as contours at -2,2,3,5 & 8 x the rms 
fluctuation in the image of it = 60;uJy . The cluster core rc = 23'.' 1 
and an area of radius 4 x rc are indicated. The positions of the 
five brightest pulsars (C, D, E, F, J) are marked with •. All other 
known pulsar positions are marked o. The size of the restoring 
beam is shown at lower left. 



of N< 30. This result is consistent with the possibility that 
there are no further radio pulsars to be detected in 47 Tuc. 

The minimum luminosity of this distribution is 
'"^"Li4()() = SrmnD^ = 0.4mJy kpc^ at 1400 MHz. Pulsar 
luminosities are often referred to 400 MHz with a spec- 
tral index of —2 assumed to allow conversion: L400 = 
^1400 (400/1400)"^. Using this approach we estimate the 
lower limit of 400 MHz luminosity of the 47 Tuc pulsar pop- 
ulation to be """L400 = 6 mJy kpc^. 



6 DISCUSSION 

In this work we have estimated two parameters of the pop- 
ulation of pulsars in 47 Tucanae: the lower bound to its 
luminosity distribution "'*"Li4oo and its size N, the total 
number of observable pulsars in the cluster. Our technique 
is similar to that used by Fruchter and Goss (1990) who 
made measurements of integrated flux densities of a num- 
ber of globular clusters and estimated the number of pulsars 
in each cluster. Fruchter and Goss used a luminosity func- 
tion corresponding to equation 3 with an assumed value of 
'"'"1/1400 ~ 0.2 mJy kpc'^, the luminosity of the weakest mil- 
lisecond pulsar known at the time of their work. In the case 
of 47 Tuc, we have used our measurement of the integrated 
flux and also the measured fluxes of the brightest members 
of the population to constrain both N and "^*"Li4oo- We 
find '"'"1/1400 ~ 0.4 mJy kpc'^, significantly greater than the 
value assumed by Fruchter and Goss. 

As we noted in the introduction, the determination of 
the radio luminosity distribution of pulsars is frustrated by 
uncertainties in pulsar distances and the difficulty of defin- 
ing a sample free of strong selection effects. To make this 
determination usefully for all currently known millisecond 
pulsars is probably impossible. The pulsars in 47 Tuc pro- 
vide the best opportunity available at present to characterise 
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0.01 0.1 1 

s (mJy) 

Figure 4. The cumulative flux distribution 5(> s) of pulsars in 
47 Tuc as measured by Camilo et al., 2000. The solid line shows 
the expected value of 5(> s) corresponding to expression 7 using 
the fitted value of the product Nsmin = 0.66mjy from Figure 
1. Corresponding values of the total pulsar number are shown on 
the right hand axis. The best estimate of total radio flux from 
the cluster core is indicated by the dashed line. The horizontal 
dotted line shows the 2ct upper bound. The vertical dotted line 
indicates the likely lower bound to the pulsar flux distribution. 



the radio luminosities of millisecond pulsars, and perhaps of 
pulsars in general. 

If wc accept the inferences made about the luminos- 
ity function of the 47 Tuc pulsar population, then the least 
luminous members of the population are nearly an order 
of magnitude brighter than the weakest pulsars in the gen- 
eral galactic population. We speculate that the origin of this 
difference relates more to the distinction between standard 
and millisecond pulsars, than to the different locations of 
the sample populations. The radio luminosity of an indi- 
vidual pulsar is expected to decline steadily, in accord with 
declining rotational energy. Ultimately the rotational energy 
reduces to the point of being insufficient to power the radio 
emission mechanism, and the pulsar crosses the "death line" 
— see for example the discussion by Chen and Ruderman 
(1993). As a class, the millisecond pulsars have very small 
values of P, so that their evolution towards the death line is 
very slow. The characteristic age at which a pulsar crosses 
the death-line (which also corresponds to a fixed mechanical 
luminosity line), can be shown to be inversely proportional 
to its magnetic field strength. Hence, compared to normal 
pulsars the age at which radio millisecond pulsars die would 
be several orders of magnitude longer. Perhaps the relatively 
high low-luminosity bound inferred for the 47 Tuc pulsars 



indicates how much additional time is required for them to 
evolve through lower luminosities to the death line. 

Most energy lost from a pulsar is radiated in the form of 
low frequency "dipole radiation" , generated by the rotation 
of the star's magnetic field and observable as the gradual in- 
crease in rotation period. A tiny fraction of the energy loss 
is visible in the radio spectrum. Real physical insight might 
be expected from successfully connecting, through measure- 
ment, the radio luminosity with the total energy loss rate. 
The pulsars in 47 Tuc present us with a measurable radio 
luminosity distribution at 1.4 GHz. Grindlay et al. (2002) 
estimated energy loss rates E of the pulsars after assuming 
that the location of each pulsar in the cluster's gravitational 
potential can be estimated from its differential dispersion 
measure. If the total radio luminosities could bo inferred 
from the 1.4 GHz measurements and the likely beam solid 
angles, some attempt could be made to connect them to the 
estimates of E. 

Finally we note that X-ray detection provides an alter- 
nate means of estimating the total population of 47 Tuc. 
Grindlay et al. (2002) have reported Chandra observations 
and place an upper limit of 90 pulsars in the cluster. The 
X-rays are attributed to thermal emission from the polar 
cap and are not beamed. Adopting the upper limits of 90 
X-ray and 30 radio pulsars, and assuming the X-rays are not 
beamed, the implied radio beaming factor is ~ 0.3. 

In conclusion, we have determined a lower bound to the 
luminosity of pulsars in the globular cluster 47 Tucanae. If 
this bound is typical of millisecond pulsars in general, then 
it would imply a considerably smaller number of pulsars in 
globular clusters than believed earlier. 
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